Equations of motion method for triplet excitation operators in graphene 
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Particle-hole continuum in Dirac sea of graphene has a unique window underneath, which in principle leaves 
a room for bound state formation in the triplet particle hole channel [Phys. Rev. Lett. 89, 016402 (2002)]. In 
this work, we construct appropriate triplet particle-hole operators, and using a repulsive Hubbard type effective 
interaction, we employ equations of motions to derive approximate eigen-value equation for such triplet op- 
erators. While the secular equation for the spin density fluctuations gives rise to an equation which is second 
order in the strength of the short range interaction, the explicit construction of the triplet operators obtained here 
shows that in terms of these operators, the second order can be factorized to two first order equations, one of 
which gives rise to a solution below the particle-hole continuum of Dirac electrons in undoped graphene. 

PACS numbers: 71.45.Gm, 81.05.ue, 81.05.uf 



INTRODUCTION 

The single particle excitations in graphene and graphite are 
characterized by a Dirac cone |nj-|4||. As for the excitations 
in the two-quasi-particle sector, adding interactions may pro- 
duce bound states in, especially particle-hole channel. Such 
bound states exhaust bosonic portion of excitation spectrum. 
In doped graphene, where an extended Fermi surface instead 
of Fermi points governs the continuum of free particle-hole 
excitations, the long range Coulomb forces binds initially free 
particle-hole pairs into spin singlet long lived bosonic exci- 
tations known as plasmons [0, 0| . Now, let us think of what 
happens in the limit where doping tends to zero? In this limit 
the area of the Fermi circle becomes smaller and smaller, so 
that the ratio of Coulomb energy to the kinetic energy in- 
creases, and the single particle picture is expected to deviate 
from the simple Dirac cone, whereby signatures of correla- 
tion effect are expected to become important in the limit of 
undoped graphene. 

The simplest model Hamiltonian which takes the domi- 
nant correlation effects into account is the Hubbard model. 
In the light of recent ab-initio estimates of the Hubbard U in 
graphene, whose unscreened value can be as large as ~ 10 
eV 10, it is important to examine possible consequences of 
such a large on-site interactions on the physical properties of 
graphene. Recently extensive quantum Monte Carlo (QMC) 
study of the phase diagram of the Hubbard model on the hon- 
eycomb lattice, suggests spin liquid ground state [0] for a 
range of U/j ~ 3 — 4, (7 being the nearest neighbor hopping 
amplitude). Therefore graphene is likely to be in the vicinity 
of a quantum spin liquid state |7[] . This scenario has been sup- 
ported by other quantum Monte Carlo studies [0] . Our recent 
QMC study suggests that the collective particle-hole excita- 
tions in sp 2 bonded planar systems are compatible with a pic- 
ture based on spin charge separation ficl . In this scenario, the 
lowest excitations are triplet states which can be interpreted 
as two-spinon bound states. It is followed by a_singlet excita- 
tion constructed from a doublon and a holon lUlll , Moreover, 
lattice gauge theory simulation of 2 + 1 dimensional QED 



predicts the critical value of the "fine structure" constant in 
graphene can be crossed in suspended graphene IU2I1 . In this 
scenario, the ground state of graphene in vacuum is expected 
to be a Mott insulator, where in the ground state, the two- 
particle sector is dominated by long-range resonating valence 
bond correlations ill 311 . Therefore, despite an intriguing sim- 
plicity of the one-particle sector of excitations in graphene, 
the two-quasi-particle sector of excitations seems to be quite 
involved and may have remarkable singlet correlations in its 
ground state. Therefore it timely to revisit the nature of spin 
excitations in undoped graphene lfl4ll from weak coupling side 
which is describe by a Dirac liquid fixed point lfl5ll . 



The collective excitation considered here, will have distinct 
features from plasmons, because: (i) Formation of plasmons 
requires doping, while here we consider undoped graphene. 
(ii) Plasmons are formed in the singlet particle-hole channel, 
as a result of long range Coulomb forces. But here we as- 
sume a short range Hubbard type interaction, and focus on the 
triplet channel of particle-hole excitations. By constructing 
equations of motion IToTl for triplet excitations formed across 
the valence and conduction band states in a Dirac cone, we 
obtain two triplet operators whose eigen-value equations are 
decoupled, and one of them displays solutions for finite val- 
ues of the short range interaction strength. We compare our 
derivation with a naive RPA-like construction of a geometric 
series lfl7ll . and show that for the triplet operators proposed 
in this work, the secular equation decouples into two first or- 
der equations in the short range interaction strength, one of 
which always does support a solution below the particle-hole 
continuum 11411 . Such a decoupling can not be achieved for 
spin density fluctuation operators IU7I1 . Since these bosonic 
excitations are not precise spin density fluctuations, their cou- 
pling to neutrons is expected to be less than the coupling of 
spin density fluctuations. We therefore discuss the coupling 
of neutrons to such excitations. 
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EFFECTIVE HAMILTONIAN 

As mentioned earlier, unlike plasmon (singlet) excitations, 
for which the long-range part of the Coulomb interaction is es- 
sential, since here we are interested in collective excitations in 
triplet (spin-flip) channel, we only need to consider the short 
range part of the interaction, as the spin-flip interactions are 
generated by short-range part of the interactions. It can be 
shown that inclusion of longer range part of the interactions 
does not lead to qualitative change in the dispersion of spin-1 
collective excitations ll 811 . Hence we start from the Hubbard 
model, 



H — Hq + Hjj 

= - 7 J2 (4a b Jc 



b )a a i<r 



Uplift 



where i, j denote sites of a honeycomb lattice, and a stands for 
spin of electrons. In this model, U ~ 10 eV is the bare value 
on-site Coulomb repulsion, and 7 ~ 2.5 eV is the hopping 
amplitude to nearest neighbor sites. To be self-contained and 
to fix the notations, we briefly summarize the change of basis 
needed to diagonalize Hq. We introduce the Fourier trans- 
forms 



where 61, 82, S3 are vectors connecting each atom in the A 
sub-lattice to its nearest neighbors. The phase of the form 
factors are defined by 7$(fc) = ^\^(k)\e lip k = e-e^s, in 
terms of which the hopping term becomes, 



*-&(u) A? Ji 

ka 



The following change of basis from (a, b) basis to (c, v) basis, 



Ha ) V2 V e 



-e- % fn 



(5) 



brings Hq to diagonal format: 



c- cr — vl vr 

ka ka ka kc 



(6) 



Operators ct and v^ a correspond to electron and hole opera- 
tors. For later reference we note the explicit relation connect- 
ing a) and operators to these basis is given by, 



1 



k 



b \ = J_ y e -ik.R j+sh l 



Ja J2 V ^a 

k 



V2 



ct 

fee 



3 ° VN^ ka ka 



(7) 



(8) 



where two atoms in the j'th unit cell are located at Rj (e A) 
and Rj+s (G B). N is the the total number of cells. The above 
Fourier expansion, transforms the non-interacting part of the 
Hamiltonian to, 



Ho = -7 E ^Haho 

with the form factor given by, 



ka 



(2) 



Now let us rewrite the short range Hubbard interaction in 
the new basis in which Hq is diagonal. The Hubbard interac- 
tion term in the exchange channel can be written as, 



Hu = -U^2 a^a^a]^ + b^b^b^b 



(9) 



$(fc) 



D ifc.(5i 



Jk.82 



Ak.53 



(3) Upon using Eqns. (0 and (O we have, 
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2N ^ 
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k,k' 



2N 



,i{p-f)-Rj + s J<Pp 



k,k' 
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(10) 

(11) 

(12) 
(13) 
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Inserting the above equations in the Hubbard term, jj J2j produces a momentum conservation constraint which can be satisfied 
by changing from k' and p 1 to a new variable q defined by 

q = k-k'=0-p, (14) 

which eventually gives 



H, 



U 

'AN 



kpq 



Expanding the Hubbard interaction in terms of electron (c) 
and hole (v) operators, generates 32 terms. Combining the 
amplitudes 1 and e lip ''^ t<p ''-i +t<Pe ^ tVi!+i from first and second 
lines of Eq. 1151 . leads to 16 types of terms with arbitrary 
number of electron and hole operators, whose amplitudes are 
of the form: 



tonian we use in this work is, 



7r~ = 

kpq 



1 ± e ivs ~ ivs -<i +iViS ~ iVi!+ ' 1 



(16) 



As can be seen from Eq. 1151 . for those terms containing 
imbalanced number of c and v operators, the amplitude of 
the process generated by interaction will be 7", while for 
those where number of conduction and valence operators are 
balanced, the interaction vertex will be proportional to 7+. 
Therefore in the long wave-length limit, |<f| — > 0, where 
7~ — » 0, we expect the following types of terms to survive 
in the effective short range interaction: 



Hi = 


u 

2N 


T- it 

^kpq Ififfl 




-qi v i>l c p+q1~ 


(17) 


H 2 = 


U 
2N 


T- it 

^kpq 'kpq 


v h c k~ 
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(18) 


H 3 = 


U 
2N 


V- if 

^kpq Ififfl 


ct cr 
fct k ~ 


-qi V li V P+^ 


(19) 


Hi = 


U 
2N 


^kpq^kpq 


v % v k- 


-qi c Pl c p+g-t 


(20) 


H 5 = 


U 
2N 


it 

^kpq 'kpq 


c % v k- 


-ql C ll V P+^ 


(21) 


H 6 = 


U 
2N 


T- it 

^kpq iftpq- 


vt cr 




(22) 



There are two more terms with their vertex strength propor- 
tional to 7 + , namely c*cc)c and v'vv'v which correspond to 
particle-hole fluctuations solely in the conduction or valence 
band, which will not contribute in the undoped graphene, as 
average occupation numbers arising from the Hartree decom- 
position of the equations of motion (see the following section) 
makes them irrelevant at this mean field level. Beyond the 
mean field, they are supposed to renormalize the bare value of 
U — > U. Therefore the tree level effective short range Hamil- 



H cS = H + Y,H a 



(23) 



The bare value of U ~ 47 is expected to get renormalized 
to a smaller value U < 2.237, beyond which an instability 
occurs 111 811 . The mean field factorization in the equation of 
motion employed here (next section) may lead to the under- 
estimation of this upper value for U, which is a known effect 
of mean field treatments lfl9ll . Therefore the physical range of 
parameters is limited to U ~ 27. More elaborate calculations 
based on exact diagonalization, as well as ab-initio quantum 
Monte Carlo calculation by us, supp orts the picture emerging 
from this effective Hamiltonian llOlj . In the following section, 
we identify appropriate triplet operators, which satisfy a sim- 
ple eigen-value equation, which correspond to singularities in 
the one-band RPA-type susceptibility. 



CONSTRUCTION OF TRIPLET OPERATORS 

As a two-band generalization of the triplet excitation in 
YBCO superconductors fl(& , consider two triplet operator de- 
fined for the particle-hole channel by, 



% V k+qV 



(24) 



These operators create triplet particle-hole excitations across 
the valence and conduction bands. Therefore, by construc- 
tion, these operators are supposed to generate (triplet) excita- 
tions in undoped graphene. To study the dynamics of these 
triplet excitations we calculate their equation of motion in a 
normal state. In the right hand side of terms generated by 
the Hubbard interaction, we perform Hartree factorization in 
terms of ap propriate occupation factors and the operator under 
study lllql . For c \^ v k+qi excitations, non-zero contributions 
are generated by, 



4 



J 

t 



fet k+qi 
A^k+qi 



u 

2N 



k+qj ^ k +qV 



2N (Sfclpr 



k,k+q,-q fet fc +<?^ 



-q V k-[ C i+li 



H + ql-^W' 



(25) 
(26) 
(27) 



where a Hartree factorization in the right hand side has been performed to generate the average occupation numbers lfl6ll . 
Similarly the non-zero contributions for triplet excitations from conduction to valence band after Hartree factorization become, 



t 

k^k+ql 



^A^k+qi 



vl.ci 



~2N [Yjklit 



k 'k,k+q-q V k^ C k+ql 



k+q) "kf k+qi 

(n 



— ft v ~ ) 

k+qi k-f' 



■2N ( Sfe ith 



ct vr , 



k 'k,k+q-q k\ k+qi 



V k+qi k\' 



(28) 
(29) 
(30) 



The above set of results can be summarized as, 



t 

k^ U k+qi 



7 f L 
k\ k+ql 



Sr. + Sr 



Si 



J c \i v k+qi 2N ^ n k+qi 



k,k+q,-q V fct fc+ ^ fct fc +9^ 



E7 



- "It) E ^5,-5 (4a~ 



I ^ V lt C k+qi 



(31) 
(32) 



Here H Q g = Y^=o ^a- Demanding right hand side of the above equations to be uiq times c| -jv k+ ^ and v ju. c k+q± < respectively, 



we obtain 



clvi 



Uq + 



k+qj °fct k+qi 



k+q) V k-\° k +qi 



U f-v 

n~ 

47V \ k 

■— (n c ~ 



— Tl~ 
k+q^ fet 



'At 



E( 1 + ^,g-C~) (4t^+w 



E^+^g-C-) (4 



fct^+g-l- 



(33) 
(34) 



With %q 



e 



oWg—Wi+s s which satisfies the property 77^" z, = 

ryS _ = 77 g _ ~ These set of equations suggest to define the 

following operators: 



0g =Efe 4«fe 



^ \ c fct^+g4 w fct c ^+9^/ ^>g 



(35) 
(36) 



where 

xo(q) -- 
xo(q) ■- 
x'M - 
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t k 




g 
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e fe+g 
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-^t 




" 

fe+gj. 
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UJq + e- k 4 
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UJ q - 




- £ fc+ 


g 







(39) 



%,g' 



%q- 



The eigen-value problem for these operators become, 



O q =U X o(q)Oq + Ux' (q)Oq, 

Oq =U X0 (q)Oq + U X 0(q)Oq, 



(37) 
(38) 



In the last equation, since 77I . 



and the fractions 



^-fc,-g' 

under the parenthesis are invariant with respect to inversion of 
the vectors k and k + q, we conclude that xo = Xo- Hence 
eigenvalue equations decouple into the following equations 
for the symmetric and antisymmetric modes constructed from 
O and O triplet operators: 



If (l - U( Xo ± Xo)) = 0. 



(40) 
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where the normalized form of our triplet operator — Oq± 
Oq is given by, 



'q 



: -u± E + v lt c Uu) 0- ± (41) 

« fe 



The normalization factor satisfies, 



(^) 2 = 4 E( 1±cos (^ 



(42) 




hui/vp 



Since we are dealing with undoped graphene, and assuming 
the temperature to be zero, the occupation numbers in the 
conduction and valence bands will be and 1, respectively. 
Therefore, the susceptibility corresponding to the above oper- 
ators reduces to, 



(43) 



Assymptotic behavior at PHC edge: 



hio — > vp\q\ 



FIG. 1: Points in fc-space located on the ellipse contribute in the 
delta function integration of S?xo. In fiu — > — ► limit the 
ellipse degenerates into two line segments where k and k' = k + q 
point in the opposite direction; 9^ — Op — ¥ n. Hence the overlap 
factors in this limit tend to a trivial factor of 1 and become irrelevant. 



In the low-energy limit where the Dirac cone linearization 
of the spectrum is valid, these integrals in the the particle-hole 
fluctuation can be analytically performed J5|,|20l]. Otherwise 
they can be computed with standard numerical procedures. 
Let us present here a geometric arguments based on a very 
peculiar constraint in the k-space which arises from the conic 
spectrum. Chiral nature of one particle eigen states, implies 
that the back-scattering will not be allowed for scattering of 
two electrons or two holes. However, for the scattering of 
a particle and a hole, the very same chiral nature according 
to which the matrix elements between a hole and an electron 
state is proportional to 1 — e i¥3c ~ lI/3£ +9, enhances the back- 
scattering between the particle and a hole. In the small if limit 
as in Fig.Q] the contributions to the imaginary part of integrals 
in Eq. d43l comes from a set of points on the ellipse defined 
by tUq = uf(|&| + \k + q\). Enhanced back-scattering in the 
particle-hole channel in the above geometry corresponds to 
the limit where ellipse degenerates into two almost parallel 
line segments, i.e. the limit uiq- ps vf|<?|. This limit corre- 
sponds to the lower edge of the particle-hole continuum where 
indeed a dominant inverse square root behavior (see Eq. l45l l 
describes the non-interacting susceptibility J5[[l3l- Such a ge- 
ometrical constraint in the phase space could be considered as 
a novel rout to 2D bosonization scheme which may find ap- 
plications in systems with Dirac cone, such as graphene and 
topological insulators 12111 . 

Now let us study the behavior of susceptibilities corre- 
sponding to the two triplet operators obtained here with the aid 
of the geometrical argument introduced above. In this limit 
the ellipse tends to two line which enforces k and k + q to be 
in opposite directions, such that, 1 ± 77^ - = l±e iV£_,¥, '=+«' — > 

1 ± e i7r . Therefore in the small q limit, Xo(o) + Xo(o) van " 
ishes, and a very large value of U will be required to excite the 



Tg~. Hence, as long as we are interested in solutions for finite 
values of U, we are left with the triplet operator 7^ whose 
eigenvalue equation can be simplified to: 



1 - UXo(q) = 0, 



(44) 



k+q , 



where we have used the symmetry of the bosonic propagator 
under the inversion symmetry (in k space) to project out the 



even part of the 1 — e 



E +9 factor. As we mentioned, 



closed form formula for the above particle-hole fluctuation 
can be obtained II5[ 12011. 



16 ^Jv A F q 2 - Ld 2 



(45) 



which despite ignoring phase factors (1 — -) is identical 
to the result obtained in Ref. 1 1-411 - The fact that ignoring 
these phase factors in Ref. 111411 in undoped graphene, does not 
change the low-energy behavior lends on the particular geom- 
etry arising from chiral nature of single-particle states (Fig.Q]). 
However, for the case of doped graphene, one has to properly 
take them into account, even for short range interactions I22I1 . 

The eigenvalue equation for is equivalent to diver- 
gence in the triplet susceptibility at random phase approxi- 
mation 123J, 



Xt R rfptt(9> W ) 



xo(q,u) 

- Uxv(q,oj)' 



(46) 
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where retarded bare susceptibility xo,aa' (q, w) in our notation 
is given by the standard particle-hole form 0], 



\ ^ V k+q 



if) (l + aa'cos(^-^)) 



2N ^ 

k 



hu — (a'ej: + ~ — aej:) + i0 + 



(47) 



Here a, a' take ± values corresponding to conduction and 
valence bands, respectively fl]. To understand the origin of 
overlap factors in this expression we note that the matrix el- 
ements of the scattering interaction V between chiral states 
(k, a) and (k', a') of the cone-like dispersion in graphene are 
given by (k',a'\V\k,a) = V(k-k') (l + aa'e^r 1 ^) /2, 
where V is the Fourier transform of the scattering potential. 
When the above phase factors are inserted into particle-hole 
bubble diagrams, give rise to the overlap factor in the free 
particle-hole propagator, Eq. (l47l i. Therefore, despite that the 
operators considered here are not exactly what one expects 
from a local spin fluctuation operator, nevertheless, the sus- 
ceptibility corresponding to them is the particle-hole bubble. 
But the main difference between spin density fluctuation and 
our triplet operators is that the former satisfies a second order 
equation in U which does not have a solution [ 17], while our 
triplet operators satisfy two first order equations, one of which 
as will be shown below has a solution below the continuum of 



free particle-hole excitations 114 1. 



SHORT RANGE VERSUS LONG RANGE INTERACTIONS 

To emphasize the importance of Eq. d40b for short range 
interactions, let us discuss how one obtains a matrix form for 
the spin density fluctuation which then leads to a second order 
secular equation lfl7ll . When the range of interactions is so 
short that the two neighboring atoms from two sub-lattices 
A, B can be resolved, an RPA like geometric series for the 
susceptibility gives rise to the following equations: 



XAA 
XAB 
XBA 
XBB 



X% + UxaaXaa + Ux%Xba 
Xab + UxaaXab + UxabXbb 
UXbaXaa + UXbbXba 
UxbaXab + UXbbXbb 



Xba 
Xbb 



where due to short range interaction U, the Hubbard inter- 
action connects x°X products in a manner that their internal 
indices are the same (Fig.|2]i. The above set of equations de- 
couples and gives two sets of determinants of the following 
form, 



I-Uxaa -Ux°a 



AB 



-Uxba l-Uxl 



BB 



0. 



(48) 



which is a quadratic equation obtained by constructing RPA- 
like series of Feynman diagrams for particle-hole pairs prop- 
agating between the lattice sites. This condition corresponds 
to the poles in the RPA susceptibility of spin density fluctua- 
tions, which will not lead to any solution (bound state) lfl7ll . 




FIG. 2: Typical geometric series for short range interactions where 
the bubbles meet either on A, or B sites. Therefore internal indices 
for short range interactions become identical. 



To see why the second order equation is peculiar to short 
range interactions, consider long range interactions, for which 
one can construct the following geometric series, 

XAA = X°AA + W ^AmXm'A 



XAB = XAB + W XAmXm'B 



XBA = X°BA + W XBmXm'A 
m,m' 

XBB = XBB + W Y XBmXm'B- 



Here due to long-range interaction W, the indices to, to' are 
not necessarily identical, and they can take both A and B val- 
ues, as the long range interaction treats both indices on the 
same footing. Summing all these equations, we can see that a 
symmetric mode decouples from the rest of equations, which 
satisfies the following equation: 



1 " W( X a AA + Xab + Xba + X%b) = 0, 



(49) 



which is first order in the (long range) interaction strength, 
W. Therefore, the factorization of the eigenvalue equation 
into first order equations in presence of short range interac- 
tions discussed for our triplet operators is not trivial, and is a 
consequence of the peculiar form of these operators. If one 
writes down the spin density operator, e.g. Eq. d551 l. one can 
see that the operators considered here do not precisely corre- 
spond to spin density fluctuations. Spin density fluctuations 
in presence of short range interactions give rise to a second 
order equation as argued above lfl7ll . However, a judicious 
choice of operators done here, manages to decouple an 
equation which otherwise is expected to be second order, into 
two first order equations. This argument not only presents the 
explicit formula for the triplet fluctuations, but also supports 
our earlier prediction of neutral triplet excitations in undoped 
graphene and graphite lfl4ll . 



NEUTRON SCATTERING CROSS SECTION 

Using Eq. d45l ) to solve the eigen-value equation ( PHI ), 
lives the following dispersion for the collective triplet exci- 
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FIG. 3: (Color online) Plot of the real part of RPA-like expression in 
undoped graphene, for U = I.97. The blank region in u — q plane 
corresponding to higher values of u> is the particle-hole continuum, 
where the imaginary part of x° becomes non-zero. As can be seen, 
a singularity in the the RPA susceptibility below the particle-hole 
continuum occurs, which indicates vanishing of the RPA denomina- 
tor, which is nothing but solution of the eigen- value equation for the 
triplet collective excitation. Energies are in units of 7. 



In polarized neutron scattering experiments, one measures, 
S{q,u) = Y J MS^\n)\ 2 S(Lo-u: n0 ) (51) 

n 

where |0) and |n) are ground and excited states of the whole 
system. As discussed in this paper, a class of approximate 
excitations are given by 



J-7T 

Af~ q 
g 



(52) 



where the normalization factor 1/Afq has been defined by 
Eq. d42l) . Contribution of this class of excitations to the struc- 
ture factor will be given by 



St r iplet(9> W ) = 



1 

If? 

g 



AT 2 
g 



(0|[Sr,7T]|0) S(u-u 9 ) 



(53) 



tations (1 



w(q) = v F q 



U 2 
32v F 



(50) 



which is valid in the limit where Dirac cone linearization ap- 
plies. When the entire band dispersion is used, one can per- 
form the integrals numerically. The solutions of the eigen- 
value equation d40b for Ta~ operator, can be visualized as 

singularities in an RPA-like susceptibility Xo/(l — Uxo)- In 
Fig. [5] we have plotted the real part of the above RPA-like ex- 
pression for few values of q and a typical value of U = 1.9j. 
The location of sharp divergences in the horizontal plane in 
this plot represents the dispersion of the neutral triplet collec- 
tive excitation generated by 7^ . 

The Dirac cone description of the electronic states of 
graphene equally holds in graphite, as long as one is interested 
in energy scales above the inter-layer hopping, t± ~ 50 meV. 
The Dirac cone description of the electronic states is valid in 
length scales much larger than the lattice spacing. Therefore, 
even in highly oriented pyrolitic graphite (HOPG) where var- 
ious planes might be slightly rotated around the z-axis with 
respect to each other, anisotropy in the momentum space can 
be safely ignored and still a Dirac cone description will remain 
valid. Hence our formulation of the spin-1 collective excita- 
tions is not only relevant to graphene, but also will be relevant 
to graphite and HOPG at energy scales above the inter-layer 
hopping. For such bulk samples, one may think of neutron 
scattering to search for the spin-1 collective excitations. How- 
ever, since the operator corresponding to the triplet excitation 
is not identical to the spin density fluctuations, the coupling of 
neutrons is expected to be renormalized by appropriate matrix 
elements. Therefore in this section, we consider the behavior 
of neutron peak intensity in the limit of small q, with explicitly 
taking our triplet operators into account. 



where we have used the fact that there are no triplet excitations 
in the ground state: (0|7^~ = 0. Moreover, note that here we 

need the vacuum expectation value of the T~ , so that the v].c^ 
term gives zero when acting on |0). Hence in the calculation 
of commutators, we drop the c\.vi part of the triplet opera- 
tor. The spin-flip operator in the present two-band situation is 
given by 



= (4i + c k) + c p-^> 



(54) 



9- = 1 v 

gf> 2-Jn 



( v k ~ c k) ( u p-«t ~ cp-^ 5 ) 



The required commutators will become 

1 V 



Sq-,b'Tg 



N 4^ 

k 

[ S q,ai 



1 



v k+qi k +gl fet fe t' 



V, 



(56) 



At zero temperature, the conduction band is empty and the 
valence band is completely filled, so that the intensity of the 
mode will be given by 



1 

AT? 



k 



q 2 



(57) 



where the asymptotic expressions for the integrals required 
above are obtained with the aid of the following expansion: 



1 



\k + q\ = k + k.q+—[q 2 - (q.k)l + 0(q 



(58) 



The oc q 2 behavior of the neutron scattering intensity makes 
the direct observation of such quanta of triplet excitations 
challenging for neutron scattering experiments. Optimum 
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spots for a neutron scattering experiments are away from the 
r point 1 1811 . Moreover, due to vanishingly small binding en- 
ergy of the triplet excitations with respect to the lower bound- 
ary of the particle-hole continuum in Fig. [3] the resulting neu- 
tron peak even if observed, maybe washed by broad spectra of 
the adjacent free particle-hole pairs. Therefore a gap opening 
mechanism in graphite (such as proximity to a superconduct- 
ing condensate, etc.) can be helpful in separating the energy 
scales associated with the expected sharp resonance peak from 
broad features associated with the continuum of incoherent 
excitations. 



SUMMARY AND DISCUSSIONS 

The secular equation obtained by Peres and coworkers jV, 
for spin density fluctuations in presence of short range inter- 
actions, is second order in U, which does not admit a solu- 
tion. However, here instead of second order equation, we ob- 
tain two set of first order equations for 7^ operators, one of 
which (7^ + ) does not lead to split-off state for finite U, while 
the other (T~~) satisfying a first order equation leads to a dis- 
persive triplet collective excitations whose energy band-width 
is on the scale of the hopping amplitude 7. Such a bosonic 
branch of excitations might be responsible for: (i) The life- 
time anomaly observed in time resolved photo-emission spec- 
troscopy of highly oriented pyrolytic graphite 112411 . (ii) The 
kink observed in the dispersion of Dirac electrons in nearly 
free standing graphene samples 12511 . (iii) The spin-flip ex- 
citations observed in artificial honeycomb lattice formed by 
quantum dots l26ll . The interpretation of such triplet mode as 
weak coupling analogue of two-spinon bound states has been 
supported by some recent Monte Carlo calculations |9|, 1 1 Qfl . 
Despite intriguing simplicity of the Dirac cone for the single- 
particle excitations of graphene/HOPG, it appears that the 
particle-hole sector of excitations is likely to be more in- 
volved, and short-range and/or Heisenberg forms of interac- 
tions maybe needed to capture the underlying singlet correla- 
tions fill. 
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